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Abstract 

Despite the abundant literature on the subject appeared in the last few years, the lattice Boltz- 
mann method (LBM) is probably the one for which a complete understanding is not yet available. 
As an example, an unsolved theoretical issue is related to the construction of a discrete kinetic 
theory which yields exactly the fluid equations, i.e., is non-asymptotic (here denoted as LB inverse 
kinetic theory). The purpose of this paper aims at investigating discrete inverse kinetic theories 
(IKT) for incompressible fluids. We intend to show that the discrete IKT can be defined in such 
a way to satisfy, in particular, the requirement of completeness, i.e., all fluid fields are expressed 
as moments of the kinetic distribution function and all hydrodynamic equations can be identified 
with suitable moment equations of an appropriate inverse kinetic equation IKE. 

PACS numbers: 47.10.ad,05.20.Dd,05.20.-y 
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I. INTRODUCTION 



Basic issues concerning the foundations classical hydrodynamics still remain unanswered. 
A remarkable aspect is related the construction of inverse kinetic theories (IKT) for hy- 
drodynamic equations in which the fluid fields are identified with suitable moments of an 
appropriate kinetic probability distribution. The topic has been the subject of theoreti- 
cal investigations both regarding the incompressible Navier-Stokes (NS) equations (INSE) 
ll [2l,j2,|j, [fa] and the quantum hydrodynamic equations associated to the Schrodinger equa- 
tion [6|. The importance of the IKT-approach for classical hydrodynamics goes beyond the 
academic interest. In fact, INSE represent a mixture of hyperbolic and elliptic pde's, which 
are extremely hard to study both analytically and numerically. As such, their investigation 
represents a challenge both for mathematical analysis and for computational fluid dynam- 
ics. The discovery of IKT l| provides, however, a new starting point for the theoretical and 
numerical investigation of INSE. In fact, an inverse kinetic theory yields, by definition, an 
exact solver for the fluid equations: all the fluid fields, including the fluid pressure p(r,t), 
are uniquely prescribed in terms of suitable momenta of the kinetic distribution function, 
solution of the kinetic equation. In the case of INSE this permits, in principle, to determine 
the evolution of the fluid fields without solving explicitly the Navier-Stokes equation, nor the 
Poisson equations for the fluid pressure [5]. Previous IKT approaches have been 

based on continuous phase-space models. However, the interesting question arises whether 
similar concepts can be adopted also to the development of discrete inverse kinetic theories 
based on the lattice Boltzmann (LB) theory. The goal of this investigation is to propose a 
novel LB theory for INSE, based on the development of an IKT with discrete velocities, here 
denoted as lattice Boltzmann inverse kinetic theory (LB-IKT). In this paper we intend to an- 
alyze the theoretical foundations and basic properties of the new approach useful to display 
its relationship with previous CFD and lattice Boltzmann methods (LBM) for incompressible 
isothermal fluids. In particular, we wish to prove that it delivers an inverse kinetic theory, 
i.e., that it realizes an exact Navier-Stokes and Poisson solver. The motivations of this work 
are related to some of the basic features of customary LB theory representing, at the same 
time, assets and weaknesses. One of the main reasons of the popularity of the LB approach 
lays in its simplicity and in the fact that it provides an approximate Poisson solver, i.e., it 
permits to advance in time the fluid fields without explicitly solving numerically the Pois- 
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son equation for the fluid pressure. However customary LB approaches can yield, at most, 
only asymptotic approximations for the fluid fields. This is because of two different reasons. 
The first one is the difficulty in the precise definition of the kinetic boundary conditions 
in customary LBM's, since sufficiently close to the boundary the form of the distribution 
function prescribed by the boundary conditions is not generally consistent with hydrody- 
namic equations. The second reason is that the kinetic description adopted implies either 



the introduction of weak compressibility 



Q, S S 3, 11, 12Lor temperature Q effects of the 



fluid or some sort of state equation for the fluid pressure [l4|]. These assumptions, although 
physically plausible, appear unacceptable from the mathematical viewpoint since they rep- 
resent a breaking of the exact fluid equations. A fundamental issue is, therefore, related to 
the construction of more accurate, or higher-order, LBM's, applicable for arbitrary values of 
the relevant physical (and asymptotic) parameters. However, the route which should permit 
to determine them is still uncertain, since the very existence of an underlying exact (and 
non-asymptotic) discrete kinetic theory, analogous to the continuous inverse kinetic theory 
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[2|, |3J], is not yet known. According to some authors [15|, H7|, H8j this should be linked to 
the discretization of the Boltzmann equation, or to the possible introduction of weakly com- 
pressible and thermal flow models. However, the first approach is not only extremely hard to 
implement 19) , since it is based on the adoption of higher-order Gauss- Hermite quadratures 
(linked to the discretization of the Boltzmann equation), but its truncations yield at most 
asymptotic theories. Other approaches, which are based on 'ad hoc' modifications of the 
fluid equations (for example, introducing compressibility and/or temperature effects 2(3]), 
by definition cannot provide exact Navier-Stokes solvers. The aim of this work is the devel- 
opment of an inverse kinetic theory for the incompressible Navier-Stokes equations (INSE) 
which, besides realizing an exact Navier-Stokes (and Poisson) solver, overcomes some of 
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the limitations of previous LBM's. Unlike Refs. [2J, |3j, where a continuous IKT was con- 
sidered, here we construct a discrete theory based on the LB velocity-space discretization. 
In such a type of approach, the kinetic description is realized by a finite number of discrete 
distribution functions fi(r, t), for % = 0, k, each associated to a prescribed discrete constant 
velocity and defined everywhere in the existence domain of the fluid fields (the open set 
Q x I ). The configuration space Q is a bounded subset of the Euclidean space R 3 and the 
time interval / is a subset of R. The kinetic theory is obtained as in 

aa 

by introducing 

an inverse kinetic equation (LB-IKE) which advances in time the distribution function and 
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by properly defining a correspondence principle, relating a set of velocity momenta with the 
relevant fluid fields. 



II. 2 - LB INVERSE KINETIC THEORY (LB-IKT) 

There are several important motivations for seeking an exact solver based on LBM. 
The lack of a theory of this type represents in fact a weak point of LB theory. Besides 
being a still unsolved theoretical issue, the problem is relevant in order to determine the 
exact relationship between the LBM's and traditional CFD schemes based on the direct 
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discretization of the Navier-Stokes equations. Following ideas recently developed [2|, y, |4J, [6( , 
we show that such a theory can be formulated by means of an inverse kinetic theory (IKT) 
with discrete velocities. By definition such an IKT should yield exactly the complete set of 
fluid equations and which, contrary to customary kinetic approaches in CFD (in particular 
LB methods), should not depend on asymptotic parameters. This implies that the inverse 
kinetic theory must also satisfy an exact closure condition. As a further condition, we require 
that the fluid equations are fulfilled independently of the initial conditions for the kinetic 
distribution function (to be properly set) and should hold for arbitrary fluid fields. The 
latter requirement is necessary since we must expect that the validity of the inverse kinetic 
theory should not be limited to a subset of possible fluid motions nor depend on special 
assumptions, like a prescribed range of Reynolds numbers. In principle a phase-space theory, 
yielding an inverse kinetic theory, may be conveniently set in terms of a quasi-probability, 
denoted as kinetic distribution function, /(x, t). A particular case of interest (investigated 
in Refs. refers to the case in which /(x, t) can actually be identified with a phase- 

space probability density. In the sequel we address both cases, showing that, to a certain 
extent, in both cases the formulation of a generic IKT can actually be treated in a similar 
fashion. This requires the introduction of an appropriate set of constitutive assumptions 
(or axioms). These concern in particular the definitions of the kinetic equation - denoted as 
inverse kinetic equation (IKE) - which advances in time /(x, t) and of the velocity momenta 
to be identified with the relevant fluid fields (correspondence principle). However, further 
assumptions, such as those involving the regularity conditions for /(x, t) and the prescription 
of its initial and boundary conditions must clearly be added. The concept [of IKT] can be 
easily extended to the case in which the kinetic distribution function takes on only discrete 
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values in velocity space. In the sequel we consider for definiteness the case of the so-called 
LB discretization, whereby - for each (r, £) £ Q x I - the kinetic distribution function is 
discrete, and in particular admits a finite set of discrete values /j(r, t) G R, for % = 0, k, each 
one corresponding to a prescribed constant discrete velocity e R 3 for % — 0, k. Let us 
now introduce the constitutive assumptions (axioms) set for the construction of a LB-IKT 
for INSE, whose form is suggested by the analogous continuous inverse kinetic theory 
The axioms, define the "generic" form of the discrete kinetic equation, its functional setting, 
the momenta of the kinetic distribution function and their initial and boundary conditions, 
are the following ones: 

1) Axiom I - LB-IKE and functional setting. Let us require that the extended fluid fields 
{V>Pi} are strong solutions of INSE, with suitable initial and boundary conditions and that 
the pseudo pressure p {t) is an arbitrary, suitably smooth, real function. In particular we 
impose that the fluid fields and the volume force belong to the minimal functional setting: 



Ve^ 3 ' 1 ^ x J), (1) 
f ie C (1 ' 0) (fi x I). 

We assume that in the set Q x I the following equation 

LD(i)fi = ^i(fi) + Si (2) 

[LB inverse kinetic equation (LB-IKE)] is satisfied identically by the discrete kinetic dis- 
tributions fi(r,t) for i = 0,k. Here fli(fi) and L D ^ are respectively the BGK and the 
differential streaming and operators, while Si is a source term to be defined. We require that 
KB-IKE is defined in the set f2 x /, so that ^li(fi) and Si are at least that C^(Cl x J) and 
continuous in Q x I. Moreover Qi(fi), to be identified as usual with the BGK operator, is 
considered for generality and will be useful for comparisons with customary LB approaches. 
We remark that the choice of the equilibrium kinetic distribution f^ q in the BGK operator 
remains completely arbitrary. We assume furthermore that in terms of fi the fluid fields 
{V,Pi} are determined by means of functional of the form M Xj [fi] = ^jfi (denoted 

i=0,8 

as discrete velocity momenta). For X = X l5 X 2 (with X\ = c 2 ,X 2 = —a*) these are re- 
lated to the fluid fields by means of the equations (correspondence principle) defined by the 
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equations Pl (r,t) - $(r) = c 2 £ fi = c 2 £ /f , and V(r,t)=£ £ a,/, = £ £ a^*, 

i=0,8 i=0,8 i=l,8 i=l,8 

where c = min{|aj| , i= 1,8} is the test particle velocity and, without loss of generality, 
fl q can be identified with a polynomial expession, with the kinetic pressure p\ replacing the 
fluid pressure p adopted previously [16|. These equations are assumed to hold identically 
in the set Q x I and by assumption, fi and f^ q belong to the same functional class of real 
functions defined so that the extended fluid fields belong to the minimal functional setting 
(Cp). Moreover, without loss of generality, we consider the D2Q9 LB discretization. 

2) Axiom II - Kinetic initial and boundary conditions. The discrete kinetic distribution 
function satisfies, for i = 0, k and for all r belonging to the closure Q, the initial conditions 
fi(r,t ) = f i(r,t ), where f i(r,t ) (for i = 0, k) is a initial distribution function defined 
in such a way to satisfy in the same set the initial conditions for the fluid fields|5i (r) = 
P {t ) + Po(j) ~ $(r) = c 2 £ f oi (r) and V G (r) = f- £ ^/^(r) To define the analogous 

i=0,8 i=l,8 

kinetic boundary conditions on 5Q, let us assume that 5Q is a smooth, possibly moving, 
surface. Let us introduce the velocity of the point of the boundary determined by the 
position vector E 5CI, defined by ~V w (r w (t),t) = 4rr w {€) and denote by n(r w ,t) the 
outward normal unit vector, orthogonal to the boundary 5Q at the point r w . Let us denote 
by f- + \r w ,t) and f^~\r w ,t) the kinetic distributions which carry the discrete velocities aj 
for which there results respectively (a, — V w ) ■ n(r w , t) > (outgoing- velocity distributions) 
and (&i — ~V W ) ■ n(r w , t) < (incoming- velocity distributions) and which are identically zero 
otherwise. We assume for definiteness that both sets, for which |a«| > 0, are non empty 
(which requires that the parameter c be suitably defined so that c > |V^|). The boundary 
conditions are obtained by suitably prescribing the incoming kinetic distribution f^~\r w , £), 
i.e., imposing (for all (r w ,t) <E 5£l x I) f^\r w ,t) = f^\r w ,t). Here f^\r w ,t) are suitable 
functions, to be assumed non- vanishing and defined only for incoming discrete velocities for 
which (a,, — V w ) ■ n(r w , t) < 0. 

3) Axiom III - Moment equations. If fi(r,t), for i = 0,k, are arbitrary solutions of 
LB- IKE [Eq.([2])] which satisfy Axioms I and II validity of Axioms I and II, we assume 
that the moment equations of the same LB-IKE, evaluated in terms of the moment op- 
erators Mx [■} = £ Xj-, with j = 1,2, coincide identically with INSE, namely that 

3 i=0,8 

there results identically [for all (r,t) e Q x I] M Xl [Ufi - tti(fi) - Si] = V ■ V = and 
M X2 [L l f l -Q l (f i )-Si} = NV = 0. 
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4) Axiom IV - Source term. The source term is required to depend on a finite number 
of momenta of the distribution function. It is assumed that these include, at most, the 
extended fluid fields {V,pi} and the kinetic tensor pressure II = 3 Y^h=o fiSaan — p D VV. 
Furthermore, we also normally require (except for the LB-IKT described in Appendix B) 
that Si(r,t) results independent of fl q (r,t), /<»(r) and f W i(r w ,t) (for i = 0,k). 

Although, the implications will made clear in the following sections, it is manifest that 
these axioms do not specify uniquely the form (and functional class) of the equilibrium 
kinetic distribution function f? q (r,t), nor of the initial and boundary kinetic distribution 
functions suitably defined. Thus, both f^ q (r,t), f i(r,t ) and the related distribution they 
still remain in principle completely arbitrary. Nevertheless, by construction, the initial and 
(Dirichlet) boundary conditions for the fluid fields are satisfied identically. In the sequel 
we show that these axioms define a (non-empty) family of parameter-dependent LB-IKT's, 
depending on two constant free parameters u c ,c > and one arbitrary real function P (t). 
The examples considered are reported respectively in the following Sec. 5,6 and in the 
Appendix B. 



III. 3 - A POSSIBLE REALIZATION: THE INTEGRAL LB-IKT 



We now show that, for arbitrary choices of the distributions fi{r,t) and f^ q (r,t) which 
fulfill axioms I- IV, an explicit (and non-unique) realization of the LB-IKT can actually be 
obtained. We prove, in particular, that a possible realization of the discrete inverse kinetic 
theory, to be denoted as integral LB-IKT, is provided by the source term [21] 



Si 



Wi 



(3) 



dpi 
dt 



a* ■ (f!-/iV 2 V- V-g + Vp) 



= Si 



where ^i^f- is denoted as first pressure term. Then the following theorem hols. 



Theorem - Integral LB-IKT 

In validity of axioms I-IV the following statements hold. For an arbitrary particular 
solution fi and for arbitrary extended fluid fields: A) if fi is a solution of LB- IKE 
[Eq.^j] the moment equations coincide identically with INSE in the set fix I; B) the initial 
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conditions and the (Dirichlet) boundary conditions for the fluid fields are satisfied identically; 
C) in validity of axiom IV the source term Si is non-uniquely defined by Eq.^j. 
Proof 

A ) We notice that by definition there results identically 

^ 1 c 2 dt [ ' 

i=0 

8 

= (5) 

i=0 

= -- [f-/iV 2 v-v ■ n + vp] 

On the other hand, by construction (Axiom I) (i — 1, fc) is defined so that there results 
identically Y^=o^i = anc ^ X^=o a *^* = 0- Hence the momenta Mx 1 ,Mx 2 of LB-IKE 
deliver respectively 

V • J2 = ( 6 ) 

i=l,8 

3^ + ^ V ' VV + V ^ + f -/^V 2 V = (7) 

i=l,8 

where the fluid fields V,pi are defined by appropriate moments [16|. Hence Eqs.flH]) and (J7J) 
coincide respectively with the isochoricity and Navier-Stokes equations. As a consequence, 
fi is a particular solution of LB-IKE iff the fluid fields {V,pi} are strong solutions of INSE. 

B) Initial and boundary conditions for the fluid fields are satisfied identically by con- 
struction thanks to Axiom II. 

C) However, even prescribing v c , c > and the real function P (t), the functional form of 
the equation cannot be unique The non uniqueness of the functional form of the source term 
5j(r, t) is assumed to be independent of f^ q (r,t) [and hence of Eq.(|2J)] is obvious. In fact, let 
us assume that Si is a particular solution for the source term which satisfies the previous 
axioms I-IV. Then, it is always possible to add to Si arbitrary terms of the form Si + 5 Si, 
with 5Si 7^ which depends only on the momenta indicated above, and gives vanishing 
contributions to the first two moment equations, namely Mx- [SSi\ = XjSSi = 0, with 

3 i=0,8 

j = 1,2. To prove the non-uniqueness of the source term Si, it is sufficient to notice that, 
for example, any term of the form 5Si = — lj F(v,t), with F(r,t) an arbitrary real 
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function (to be assumed, thanks to Axiom IV, a linear function of the fluid velocity), gives 
vanishing contributions to the momenta Mx 1 , Mx 2 ■ Hence Si is non-unique. 

The implications of the theorem are straightforward. First, manifestly, it holds also in 
the case in which the BGK operator vanishes identically. This occurs letting v c = in the 
whole domain Q x I. Hence the inverse kinetic equation holds independently of the specific 
definition of /f 9 (r,i). 

An interesting feature of the present approach lies in the choice of the boundary condition 
adopted for /j(r,i), which is different from that usually adopted in LBM's [see for example 
(qI for a review on the subject]. In particular, the choice adopted is the simplest permitting 
to fulfill the Dirichlet boundary conditions [imposed on the fluid fields]. This is obtained 
prescribing the functional form of fi(r,t) on the boundary of the fluid domain (5Q), which 
is identified with a function f i(r, t). 

Second, the functional class of /j(r,t), f^ g (r,t) and of f i(r, t) remains essentially arbitrary. 
Thus, in particular, the initial and boundary conditions, specified by the same function 
/oi(r, t), can be suitably defined. As further basic consequence, f^ q (r,t) and fi(r,t) need 
not necessarily be Galilei-invariant (in particular they may not be invariant with respect 
to velocity translations), although the fluid equations must be necessarily fully Galilei- 
covariant. As a consequence it is always possible to select f^ q (r,t) and f D i(r,t) based on 
convenience and mathematical simplicity. Thus, besides distributions which are Galilei 
invariant and satisfy a principle of maximum entropy (see for example 22|, [23]), it is always 
possible to identify them [i.e., f^ q (r,t),f i(r,t)] with a non-Galilean invariant polynomial 
distribution. We mention that the non-uniqueness of the source term S% can be exploited 
also by imposing that / 4 e<? (r,£) results a particular solution of the inverse kinetic equation 
Eq.(j2J) and there results also f oi (r,t) = f^ q {r,t). 



IV. 4 - CONCLUSIONS 



In this paper we have presented the theoretical foundations of a new phase-space model for 
incompressible isothermal fluids, based on a generalization of customary lattice Boltzmann 
approaches. We have shown that many of the limitations of traditional (asymptotic) LBM's 
can be overcome. As a main result, we have proven that the LB-IKT can be developed in 
such a way that it furnishes exact Navier-Stokes and Poisson solvers, i.e., it is - in a proper 
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sense - an inverse kinetic theory for INSE. The theory exhibits several features, in particular 
we have proven that the integral LB-IKT (see Sec. 3): 

1. determines uniquely the fluid pressure p(r, t) via the discrete kinetic distribution func- 
tion without solving explicitly (i.e., numerically) the Poisson equation for the fluid 
pressure. Although analogous to traditional LBM's, this is interesting since it is 
achieved without introducing compressibility and/or thermal effects. In particular 
the present theory does not rely on a state equation for the fluid pressure. 

2. is complete, namely all fluid fields are expressed as momenta of the distribution func- 
tion and all hydrodynamic equations are identified with suitable moment equations of 
the LB inverse kinetic equation. 

3. allows arbitrary initial and boundary conditions for the fluid fields. 

4. is self- consistent: the kinetic theory holds for arbitrary, suitably smooth initial condi- 
tions for the kinetic distribution function. In other words, the initial kinetic distribu- 
tion function must remain arbitrary even if a suitable set of its momenta are prescribed 
at the initial time. 

5. the associated the kinetic and equilibrium distribution functions can always be chosen 
to belong to the class of non-Galilei-invariant distributions. In particular the equilib- 
rium kinetic distribution can always be identified with a polynomial of second degree 
in the velocity. 

6. is non- asymptotic, i.e., unlike traditional LBM's it does not depend on any small 
parameter, in particular it holds for finite Mach numbers. 

The main result of the paper is represented by the construction of an explicit realization of 
the LB-IKT for the incompressible Navier-Stokes equations. The construction of a discrete 
inverse kinetic theory of this type for the incompressible Navier-Stokes equations represents 
an exciting development for the phase-space description of fluid dynamics, providing a new 
starting point for theoretical and numerical investigations based on LB theory. In our view, 
the route to more accurate, higher-order LBM's, here pointed out, will be important in order 
to achieve substantial improvements in the efficiency of LBM's in the near future. 
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